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ABSTRACT 

We have fitted the Chandra Low Energy Transmission Grating spectrum of SS Cygni in outburst with a single 
temperature blackbody suffering the photoelectric opacity of a neutral column density and the scattering opacity of 
an outflowing wind. We find that this simple model is capable of reproducing the essential features of the observed 



spectrum with the blackbody temperature T\,\ « 250 ± 50 kK, hydrogen column density Na « 5.0lj ? x 10 ly cm 



119, 



fractional emitting area / m 5.6^ 5 x 10 , boundary layer luminosity Lu ~ 5^ x 10 erg s , wind velocity 
V « 2500 km s" 1 , wind mass-loss rate M w w 1.1 X 10 16 g s -1 , and arbitrary values of the wind ionization fractions 
of 20 ions of O, Ne, Mg, Si, S, and Fe. Given that in outburst the accretion disk luminosity Ldi.sk ~ 1 x 10 35 erg s" 1 , 
Ly/Ldisk ~ O.O5+Q03, which can be explained if the white dwarf (or an equatorial belt thereon) is rotating with an 



angular velocity f2 w< j w 0.7^ 2 Hz, hence V rot sin i ~ 2300 km s 

Subject headings: accretion, accretion disks — binaries: close - 
individual (SS Cygni) — X-rays: binaries 



novae, cataclysmic variables — stars: 



1. INTRODUCTION 

According to simple theory, the boundary layer between the 
accretion disk and the surface of the white dwarf is the dom- 
inant source of high energy radiation in nonmagnetic cata- 
clysmic variables (CVs). Unless the white dwarf is rotating 
rapidly, the boundary layer luminosity L b \ w GM W( jM/2/? w( j, 
where M w d and /? wc | are respectively the mass and radius of 
the white dwarf and M is the mass-accretion rate. When M is 
low, as in dwarf nova in quiescence, the boundary layer is op- 
tically thin to its own radiation and its temperature is of order 
the virial temperature T wu = GM^m-a/^kR^i ~ 10 keV. When 
M is high, as in nova-like variables and dwarf novae in out- 
burst, the boundary layer is optically thick to its own radia- 
tion and its temperature is of order the blackbody temperature 



> ■ T bb = (GM^M/SnafRl^ 4 ~ 100 kK. 



tra of such systems should be dominated by broad lines. 

Consistent with this picture, in outburst the EUV spectra of 
the high inclination dwarf novae OY Car (i « 83°; Mauche & 
Raymond 2000) and WZ Sge (i « 75°; Wheatley & Mauche 
2004, in preparation) are dominated by broad lines, while that 
of the low inclination dwarf novae VW Hyi (i ~ 60°; Mauche 
1996) and SS Cyg (i - 40°; Mauche, Raymond, & Mattei 1995) 
are dominated by the continuum. U Gem (i w 70°; Long et al. 
1996) is a transitional case, as its EUV flux is partially eclipsed 
and its EUV spectrum contains of broad emission lines super- 
posed on a bright continuum. 

The general case of the radiation transfer of the boundary 
layer flux through the accretion disk wind is complicated, but 
schematically we expect that in high inclination systems the 
EUV spectrum will be of the form F\e~ Tx (1 -e~ Tx ), while in 
low inclination systems it will be of the form Fx e~ Tx e~ Tx , where 
Fx is the intrinsic boundary layer spectrum, t' x is the photoelec- 
tric optical depth of the wind and the interstellar medium, and 
tx is the scattering optical depth of the wind. Mauche & Ray- 
mond (2000) showed that the second version of this model de- 
scribes the main features of the EUVE spectrum of OY Car in 
superoutburst, and we show here that the first version of this 
model describes the main features of the Chandra Low En- 
ergy Transmission Grating (LETG) spectrum of SS Cyg in out- 
burst. A preliminary discussion of this spectrum is supplied 
by Mauche (2004), while the quasi-coherent oscillations of the 
EUV flux is discussed by Mauche (2002). 

2. CHANDRA OBSERVATION AND DATA ANALYSIS 

Our pre-approved target-of-opportunity Chandra LETG plus 
High Resolution Camera (HRC) observation of SS Cyg was 
performed between 2001 lanuary 16 21 h 13 m and January 17 
10 h 50 m UT during the plateau phase of a wide normal (asym- 
metric) dwarf nova outburst that began on January 12 (Fig. 1). 
When the observation began, SS Cyg had been at maximum op- 
tical light (V w 8.5) for approximately 3^ days; long enough, 

1 This paper is dedicated to the memory and accomplishments of my colleague and friend Janet Akyiiz Mattei, who died on 2004 March 22 after a long battle with 
acute myelogenous leukemia. Her passing is a great loss to the astronomical community, both amateur and professional. 



The compact nature of the boundary layer is established by 
the narrow eclipses of the hard X-ray flux of high inclination 
dwarf novae in quiescence (Wood et al. 1995; Mukai et al. 
1997; van Teeseling 1997; Pratt et al. 1999a; Wheatley & West 
2003) and the short periods [P ~ 10 s ~ 2 7 r(/? 3 vd /GM wd ) 1 / 2 ] of 
the oscillations of the soft X-ray and extreme ultraviolet (EUV) 
flux of nova-like variables and dwarf novae in outburst (Warner 
2004; Mauche 2004). However, an extended source of soft X- 
ray and EUV flux is required to explain the lack of eclipses in 
high inclination nova-like variables and dwarf novae in outburst 
(Naylor et al. 1988; Wood, Naylor, & Marsh 1995; Pratt et al. 
1999b; Mauche & Raymond 2000; Pratt et al. 2004). 

Using the Extreme Ultraviolet Explorer {EUVE) light curve 
and spectrum of OY Car in superoutburst, Mauche & Raymond 
(2000) argued that the source of the extended EUV flux of high- 
M nonmagnetic CVs is the accretion disk wind. In this picture, 
in high inclination systems the boundary layer is hidden by the 
accretion disk for all binary phases, but its flux is scattered into 
the line of sight by resonance transitions of ions in the wind. 
Because the scattering optical depths of these transitions are 
much higher than the Thompson optical depth, the EUV spec- 
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given the delay of approximately 1 i days between the rise of 
the optical and EUV light curves of SS Cyg (Mauche, Mattei, 
& Bateson 2001), for the EUV flux to have reached maximum 
and for the system to have reached a quasi-steady state. In- 
deed, the zero- and ± first-order LETG/HRC count rates were 
roughly constant throughout the 47.3 ks observation at 1.5 and 
4.9 counts s" 1 , respectively. 
7 
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FIG. 1.— AAVSO light curve of SS Cyg (Mattei 2003). Black rectangle 
marks the interval of the Chandra LETG observation. 

The files used for analysis were created on 2002 January 7 by 
the pipeline data reduction software using CI AO 2.0 with AS- 
DCS version 6.5.1. Events were extracted from the level 2 evt 
file, and the mean spectrum and source and background region 
masks were extracted from the level 2 pha file. The observation- 
average fluxed spectrum (units of erg cm" 2 s" 1 A" 1 ) was derived 
from the ± first-order pha source and background count spectra 
and the "canned" ± first-order LETG/HRC effective area files 
distributed with Chandra CALBD 2.21. To avoid the HRC chip 
gaps, the wave bands A = 60.2-67.2 A for the positive spectral 
order and A = 49.9-57.2 A for the negative spectral order were 
excluded. The spectrum used for subsequent analysis combines 
± first orders and is binned up in wavelength by a factor of four 
to AA = 0.05 A. 

As shown by Mauche (2004), the resulting LETG spectrum 
contains two distinct components: (1) an X-ray component 
shortward of 42 A consisting of a bremsstrahlung continuum 
and emission lines of H- and He-like C, N, O, Ne, Mg, and 
Si and L-shell Fe (most prominently Fe XVII) and (2) an EUV 
component extending from 42 A to 130 A that consists of a 
continuum significantly modified by a forest of broad absorp- 
tion features (see Fig. 3). In the 72-130 A wave band, the 
LETG spectrum is essentially identical to that measured by the 
EUVE short wavelength spectrometer during the 1993 August 
(Mauche, Raymond, & Mattei 1995) and 1996 October (Wheat- 
ley, Mauche, & Mattei 2003) outbursts of SS Cyg. The X-ray 
component of the LETG spectrum is not discussed here further, 
but is understood to be due to the residual optically thin portion 
of the boundary layer. 

3. SPECTRAL ANALYSIS 

To fit the EUV spectrum of SS Cyg, we used the simple 
model described by Mauche & Raymond (2000), which as- 
sumes that the boundary layer spectrum is that of a single 
temperature blackbody suffering the photoelectric opacity of 
a neutral column density and the scattering opacity of an out- 
flowing wind. The boundary layer spectrum is then specified 
by the blackbody temperature 7m, fractional emitting area / 
(hence luminosity Lu = 47r//? 2 vd cr7y), and hydrogen column 



density Nh (assuming the Rumph, Bowyer, & Vennes 1994 
cross section with H I, He I, and He II with abundance ratios 
of 1:0.1:0.01). The properties of the wind are specified by its 
mass-loss rate Mo = 47rr 2 ^m H nv = 7.0 x 10 15 g s -1 via fiducial 
values of the wind radius r = 10 10 cm, density n = 10 10 cm -3 , 
and velocity v = 2500 km s" 1 , where the value of the wind ve- 
locity was chosen to reproduce the observed blueshift of the 
absorption lines. In the Sobolev approximation, the wind op- 
tical depth t\ = (ne 2 /m s c)nAXXijfij (gi/Y>g\)(dv/ dr)~ l , where 
dv/dr = 2500 km s" 1 /10 10 cm = 0.025 s" 1 is the wind velocity 
gradient, Ay, /ij, and gi are respectively the wavelengths, oscil- 
lator strengths, and statistical weights of the various resonance 
transitions (Verner, Verner, & Ferland 1996), A are the elemen- 
tal abundances relative to hydrogen (Anders & Grevesse 1989), 
X are the ionization fractions, and the other symbols have their 
usual meanings. To account for the velocity of the wind along 
the line of sight, the transition wavelengths were blueshifted by 
AAy = Aijv/c. To match the observed widths of the absorption 
lines, the opacity for all lines was distributed as a Gaussian with 
FWHM =1.0 A. Finally, we assumed that the source distance 
d = 160 pc (Harrison et al. 2000) and the white dwarf mass 
M„d = 1.0 M B , hence the white dwarf radius 7\ wc j = 5.5 x 10 8 



cm. 



With these assumptions, the free parameters of the model 
are the boundary layer temperature 7m, fractional emitting area 
/, hydrogen column density Nh, wind mass-loss rate M w = 
m x Mo, and ionization fractions X. Operationally, the model 
was fit to the data between 45 A and 125 A using a grid search 
in the variables 7m, Nh, and m, with X = 1 and / determined 
by the model normalization that minimizes x 2 between the 
model and the data. The resulting best-fit model has x 2 P er 
degree of freedom (dof) = 5961/1596 = 3.74 with T bi = 250 
kK, N H = 5.0 x 10 19 cm" 2 , / = 5.0 x 10"\ and m = 1.6 and 
contains 34 lines of O V-O VI, Ne V-Ne VIII, Mg V-Mg X, 
Si V-Si XI, S VII-S X, Fe VII-Fe X, and Fe XII-Fe XV (3 1 ions 
of six elements) that depress the continuum by at least 10%. 
Next, 7m, Nn, and m were frozen and X was allowed to vary 
for each of these 31 ions. Table 1 shows the evolution of x 2 as 
sequential steps were taken in this new parameter space. The 
resulting best-fit model has x 2 /dof = 3333/1565 = 2.13 with 
/ = 5.6 x 10" 3 and 41 lines of O V-O VI, Ne V-Ne VIII, Mg V- 
Mg IX, Si Vl-Si IX, S VII, Fe VII, Fe IX, Fe XIII, and Fe XV 
(20 ions of six elements) that depress the continuum by at least 
10%. The resulting contours of x 2 and / in the [7m, Nh] param- 
eter plane are shown in Figure 2 and the best-fit model is shown 
superposed on the data in Figure 3. 

With x 2 /dof = 2.13, the fit shown in Figure 3 is unacceptable 
in a statistical sense, but it clearly reproduces the essential fea- 
tures of the observed spectrum. Ignoring the unmodeled C V 
Hea emission feature at 41 A, the largest systematic residuals 
between the model and the data are in the neighborhood of 48, 
65, 86, and 92 A, where the data is above the model, but this 
discrepancy probably could be accommodated in a more de- 
tailed model by raising the continuum and increasing the wind 
optical depth. The goodness of the fit is due in large part to the 
arbitrary variations of the ionization fractions, which decrease 
X 2 /dof from 5961/1596 = 3.74 to 3333/1565 = 2.13 after 14 
steps in x 2 space. The reasons for the sequential improvements 
in the fit documented in Table 1 are often clear [e.g., increasing 
X(Mg VIII) increases the depth of the model absorption fea- 
tures at w 69.0 A and « 74.4 A; increasing X(Si VII) increases 
the depth of the model absorption features at w 67.6 A and 
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« 69.0 A; zeroing outX(Si V) removes a strong model absorp- 
tion feature at w 96.4 A], but others are less so [e.g., increasing 
X(Fe VII) increases the depth of a number of absorption fea- 
tures longward of « 102 A]. While the individual values of the 
ionization fractions in Table 1 are not unreasonable, as a set 
they make little sense physically. 



TABLE 1 
Model Ionization Fractions 



Ion 


X 


x 1 


Mg VIII 


4.0 


5210 


Si VII 


3.2 


4742 


FeVII 


10. 


4276 


Si V 


0.0 


4089 


MgVII 


1.6 


3938 


FeXIII 


4.0 


3775 


Mg V 


2.5 


3660 


Mg VI 


0.5 


3567 


FeXII 


0.0 


3511 


S IX 


0.0 


3464 


Ne VIII 


0.6 


3418 


Si X 


0.0 


3383 


Mg VIII 


3.2 


3352 


Si VI 


1.3 


3333 
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FIG. 2. — Contours of \ 2 = [l-5.3.6]x^ in (gray curves) and log/ (black 
curves) in the [TbhMt] parameter plane for the model described in §3. Re- 
gions of parameter space shaded gray are excluded by the requirements that 
Nft > 3.5 X 10 19 cnT 2 and / < 1. Cross marks the parameters of the model 
shown in Fig. 3. 

To see this, we used XSTAR vl.07 to calculated the ioniza- 
tion fractions of a photoionized plasma for various values of 
the ionization parameter £ = Lu/nr 2 . With Lbi = 47r/7? 2 vd <77^ = 
4.7 x 10 33 erg s" 1 and the fiducial values of n and r, log£ = 3.7. 
Assuming that the plasma is optically thin and that the pho- 
toionizing source is a 250 kK blackbody, the XSTAR runs 
demonstrate that the ionization parameter must span the range 
-1 < log£ < 5 to account for the range of ions present in the 
model. 

Given this result, we investigated whether it is possible to 
obtain good fits to the data with less arbitrary sets of ionization 
fractions. In the first attempt, we fixed the model parameters 
7bi and Nh at their nominal values, set the ionization fractions 
equal to their peak values for -1 < log£ < 5, and adjusted m 
and /. This model produced only a small improvement in the 
fit relative to that with X = 1, since the peak ionization frac- 
tions differ by factors of only a few. In the second attempt, we 
weighted the ionization fraction distributions X(£) by the func- 
tion e ~^ n ' in ^e~^^" r ", which suppresses the ionization fractions 



at the lower and upper ends of the range. However, varying the 
exponential cutoffs £ m ; n and £ max and m and / again produced 
only a small improvement in the fit relative to that with X = 1 . 
We conclude that our simple model is capable of producing rea- 
sonably good fits to the data only if the ionization fractions are 
allowed to vary arbitrarily. 

4. DISCUSSION 

We have fitted the Chandra LETG spectrum of SS Cyg in 
outburst with a single temperature blackbody suffering the pho- 
toelectric opacity of a neutral column density and the scattering 
opacity of an outflowing wind. The best-fit model, shown in 
Figure 3, has parameters T hx = 250 kK, N H = 5.0 x 10 19 cm -2 , 
f = 5.6 x 10" 3 , M w = 1.1 x 10 16 g s" 1 , and the values of X listed 
in Table 1 (X = 1 otherwise). Assuming that the wind is opti- 
cally thin and photoionized by a 250 kK blackbody, the range of 
ions present in the model requires that the ionization parameter 
-1 < log£ < 5. 

Given these results, in outburst the boundary layer luminos- 
ity Lbi = 4.7 x 10 33 erg s" 1 . A lower limit on the value of this 
quantity is imposed by the value of the interstellar H column 
density, which Mauche, Raymond, & Cordova (1988) estimated 
to be 3.5 x 10 19 cm" 2 based on the curve-of-growth of interstel- 
lar absorption lines in high resolution International Ultraviolet 
Explorer spectra of SS Cyg in outburst. For -1 < log£ ^ 5, 
H and He are fully ionized, so the wind should have no photo- 
electric opacity in the EUV bandpass, and Nh should be equal 
to the interstellar H column density. With N H = 3.5 x 10 19 cm" 2 
and 7bi = 300 kK, Figure 2 shows that / = 1.1 x 10" 3 , hence 
Lbi = 1.7 x 10 33 erg s" 1 . At the other end of the x 2 ellipse, 
with Tbi = 200 kK and N H = 7.9 x 10 19 cm" 2 , / = 6.6 x 10" 2 , 
hence Lbi = 2.3 x 10 34 erg s" 1 . We conclude that a reason- 
able estimate of the boundary layer luminosity Lbi ~ 5^ 8 x 
10 33 (J/160 pc) 2 ergs" 1 . 

How does this compare with the accretion disk luminosity? 
Polidan & Holberg (1984) used optical through far ultravio- 
let spectra of SS Cyg in outburst to determine that L d i s k « 
1 x 10 35 (d/160 pc) 2 erg s" 1 , so L bl /L disk w 0.05^. Theoret- 
ically, we expect that this ratio is equal to one unless the white 
dwarf is rotating rapidly, in which case Lbi/Ldi s k = (1-lu) 2 , 
where ui = fi W d/^K(^wd)> ^wd = 27r/f S pin is the angular veloc- 
ity of the white dwarf, and ft K (fl w d) = {GM^/Rl^l 2 w 0.9 
Hz is the Keplerian angular velocity of material just above its 
surface (Popham & Narayan 1995). The white dwarf can be 
spun up to the required rate by accreting an amount of mass 
AM = M w d-M w d.i, where M w d.i is the initial white dwarf mass, 
M w d = M„d,i (1 -4fe 2 w/3)" 3 ' 4 is the final white dwarf mass, and 
k is the radius of gyration (Langer et al. 2000). With k w 0.4 and 
uj = 1 -(L bl /L dis k) 1/2 » 0.78^, AM « 0.13M wd - 0.1 M Q . If 
rapid rotation is the cause of SS Cyg's low boundary layer lu- 
minosity, the white dwarf angular velocity f2 w d = cuflKiRwd) ~ 
0.7±Q2 Hz, hence the spin period P sp ; n w 9^ 4 s (comparable 
to the period of the quasi-coherent oscillations), the rotation 
velocity V rot = ri w d^wd ~ 3800^ 4 3q km s" 1 , and V rot sin/ ~ 
2300 km s" 1 . This last quantity is rather large, given that 
V rot sin ; w 1200 km s" 1 for WZ Sge and V rot sin i < 400 km s" 1 
for six other nonmagnetic CVs (Sion 1999). 

Finally, we note that in the single scattering limit, conser- 
vation of momentum limits the mass-loss rate of a radiatively 
driven wind to M max = L/VoqC, where L is the system luminos- 
ity and Voc is the wind terminal velocity. For SS Cyg in out- 
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FIG. 3. — Chandra LETG/HRC spectrum and best-tit model spectrum of SS Cyg in outburst. Data are shown by the black histogram, the 1 <r error vector by the 
blue histogram, the absorbed blackbody continuum by the smooth blue curve, and the net model spectrum by the thick red curve. In panels (b) and (c) the individual 
ion spectra are shown by the blue curves and the strongest lines in the model are labeled. 
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burst, L = L disk +L b i ~ 1 x 10 35 erg s" 1 and Voo w 3000 km s" 1 , 
so M max « 1 x 10 g s -1 . We derived a wind mass-loss rate 
M w w 1 x 10 16 g s" 1 for SS Cyg in outburst, so M w /M max w 1. 
In r Sco (BO V), P Cyg (Bl Ia + ), C Pup (04 f), and e Ori (BO 
la) M w /M max = 0.02, 0.22, 0.33, and 0.65, but in WR1 (WN5) 
M w /M m - ix = 60 (Lamers & Cassinelli 1999), so we cannot con- 
clude from this result that a mechanism other than radiation 
pressure is needed to drive the wind of SS Cyg. 

5. SUMMARY 

We have fitted the Chandra LETG spectrum of SS Cyg in 
outburst with a single temperature blackbody suffering the pho- 
toelectric opacity of a neutral column density and the scat- 
tering opacity of an outflowing wind. Figures 2 and 3 show 
that this simple model is capable of reproducing the essen- 
tial features of the observed spectrum with 7bi ~ 250 ± 50 kK, 
N H m 5.(tf;| x 10 19 cm" 2 , / « 5.6!f 5 x 10" 3 , L bl « 5^ 8 x 
10 33 erg s" 1 , and the wind velocity v w 2500 km s" 1 , mass-loss 
rate M w w 1.1 x 10 16 g s" 1 , and ionization fractions X listed in 
Table 1 (X = 1 otherwise). Given that in outburst the accretion 
disk luminosity L d isk ~ 1 x 10 35 erg s" 1 , L b i/L disiL w O.OS^q®, 
which can be explained if the white dwarf (or an equatorial belt 
thereon) is rotating with an angular velocity f2 w d ~ 0.7^ 2 Hz, 
hence V ro t sin ~ 2300 km s" 1 . 

We have now used our simple one dimensional spectral 
model to fit the EUVE spectrum of OY Car in superoutburst 
(Mauche & Raymond 2000) and the Chandra LETG spectra 
of WZ Sge in superoutburst (Wheatley & Mauche 2004, in 
preparation) and SS Cyg in outburst. Although the model is 
unable to reproduce the asymmetries of the EUV emission fea- 
tures of OY Car and WZ Sge, or to produce an acceptable fit 



to the EUV spectrum of SS Cyg even with arbitrary values of 
the ionization fractions, we believe that it captures the essential 
physics responsible for shaping the EUV spectra of these sys- 
tems, and that it returns reasonable estimates of the values of 
the physical parameters that describe the boundary layer (effec- 
tive temperature, emitting area, luminosity) and wind (velocity 
and mass-loss rate) of high-M CVs. To make further progress in 
this area, we need to expand the model from one dimension to 
three dimensions; account for the velocity, density, and ioniza- 
tion structure of the wind; and compute the paths of boundary 
layer photons as they scatter in the wind. This could be accom- 
plished with a Monte Carlo radiation transfer code, which we 
are in the process of developing (Mauche et al. 2004). 

Our pre-approved Chandra target-of-opportunity observa- 
tion of SS Cyg was made possible by the optical monitoring 
and alerts provided by the members, staff (particularly E. Waa- 
gen), and director, J. Mattei, of the AAVSO, and by the ef- 
forts of Chandra X-Ray Observatory Center Director H. Tanan- 
baum, Mission Planner K. Delain, and the Chandra Flight Op- 
erations Team at MIT. We acknowledge with thanks the vari- 
able star observations from the AAVSO International Database 
contributed by observers worldwide and used in this research. 
Support for this work was provided by NASA through Chan- 
dra Award Number GO 1-2023 A issued by the Chandra X-Ray 
Observatory Center, which is operated by the Smithsonian As- 
trophysical Observatory for and on behalf of NASA under con- 
tract NAS8-39073. This work was performed under the aus- 
pices of the US Department of Energy by University of Califor- 
nia, Lawrence Livermore National Laboratory under contract 
W-7405-Eng-48. 
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